Nanoplasmonic Lattices for Ultracold Atoms 
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We propose to use sub- wavelength confinement ol light associated with the near field of plasmonic 
systems to create nanoscale optical lattices for ultracold atoms. Our approach combines the unique 
coherence properties of isolated atoms with the sub- wavelength manipulation and strong light- matter 
interaction associated with nano-plasmonic systems. It allows one to considerably increase the 
energy scales in the realization of Hubbard models and to engineer effective long-range interactions 
in coherent and dissipative many-body dynamics. Realistic imperfections and potential applications 
are discussed. 
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Coherent optical fields provide a powerful tool for ma- 
nipulating ultracold atoms [TJ [2J. However, diffraction 
sets a fundamental limit for the length-scale of such ma- 
nipulations, given by the wavelength of light [3] . In par- 
ticular, the large period of optical lattices determines the 
energy scale of the associated many-body atomic states 
[3HS]- The resulting scaling can be best understood by 
noting that in the first Bloch band the maximum atomic 
momentum ~ l/£, where I is the lattice spacing. This 
sets the maximum kinetic energy to h 2 /m£ 2 [7j. For 
conventional optical lattices the lattice spacing is set by 
half the wavelength of the trapping light ~ 500 nm; this 
yields corresponding tunneling rates of up to a few tens 
of kHz. Additionally, for atoms in their electronic ground 
states interactions are restricted to short range. 

In this Letter we propose and analyze a novel approach 
to the realization of high-density optical lattices using 
the optical potential formed from the near field scat- 
tering of light by an array of plasmonic nanoparticles. 
By bringing atom trapping into the subwavelength and 
nanoscale regime we show that the intrinsic scales of tun- 
neling and onsite interaction for the Hubbard model can 
be increased by several orders of magnitude compared to 
conventional optical lattices. In addition, subwavelength 
confinement of the atoms results in strong radiative inter- 
actions with the plasmonic modes of the nanoparticles [5] . 
The coupled atom-plasmon system can be considered as 
a scalable cavity array that results in strong, long range 
spin-spin interactions between the atoms with both dissi- 
pative and coherent contributions HH] • Such a system 
can be used for entanglement of remote atoms as well as 
for novel realizations of coherent and dissipative many- 
body systems. 

To illustrate our approach we first consider a single 
metallic nanosphere in vacuum illuminated by a plane 
wave. For spheres small compared to a wavelength the 
dominant contribution to the scattered field is the dipole 



term, where the induced dipole moment in the sphere is 
given by p = a(ui)Eo with 
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where a is the radius of the sphere and e is the permitivity 
of the nanosphere '11.. The total electric field is 
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Near e(ui sp ) = —2 there is a plasmon resonance and 
the scattered field can be engineered to create an op- 
tical dipole trap as depicted in Fig. lb. Specifically, 
when the applied field is linearly polarized on the blue 
side of the plasmon resonance then the induced dipole 
will be ~ 7T out phase with the incident field, leading 
to two intensity minima along the polarization direc- 
tion at the positions Zj, = ±2a 3 uj 2 p /(uj 2 — ui 2 p ), where 
we took a Lorentzian polarizability near the resonance 
a(ui) — 4ireQa 3 uj 2 p /(uj 2 p — uj 2 — iu!K), with k the linewidth. 
Similarly, for red detuned, circularly polarized light, 
there will be two minima along the propagation axis. An 
atom can be trapped in these intensity minima via op- 
tical dipole forces provided the field is also on the blue 
side of the atomic resonance [2J. The trapping potential 
is given by Ml 2 /S, where O = /xq • E/H is the Rabi fre- 
quency, fio is the atomic dipole moment, and S — u a — uj 
is the detuning between the atom and laser. Expand- 
ing near the trap minima gives the trapping frequency 

w t = 9^Re(a) 2 /M 2 ~ hn 2 /Sma 2 . 

The trap depth can be controlled by applying a second 
field with the opposite polarization because this induced 
dipole constructively interferes with the applied field at 
the position of the trap, reducing the overall trap depth 
as illustrated in Fig. lc. Using this method, the atoms 
can be loaded into the near field traps by starting with 
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FIG. 1: a) Illustration of the relevant physics in the plasmonic lattice, b) Illustration of how to engineer a blue-detuned optical 
dipole trap by driving on the blue side of the plasmon resonance, c) Atomic potential for Rb including van der Waals (vdw) 
for trapping above a single silver nanoshell. Dotted line shows how to weaken the trap by applying circularly polarized light 
perpendicular to the trapping light. (Inset) Real (dashed) and imaginary (solid) part of the dipole polarizability for a sphere 
and the nanoshell with a 15 nm radius and 13.85 nm Si02 core, d) y-z contours of atomic potential in MHz for a line of nine 
spheres in the center of a 45x45 square lattice with a 60 nm lattice spacing, black regions are where the potential is negative 
due to vdw, spheres are shown in white. The nanoshells are silver with a 15 nm radius and 13.65 nm Si02 core, the trapping 
light is red detuned (704 nm) wrt to the plasmon resonance (682 nm) and applied from above with rotating x-y polarized light. 



a cold, dense gas of atoms in a large trap and then adia- 
batically turning on the near field traps. 

We now address several practical considerations. First 
of all, for alkali atoms there is a large disparity between 
the natural plasmon resonance and the atomic trapping 
transitions. For a solid silver sphere the plasmon reso- 
nance occurs near 350 nm |12) . compared to 780 nm for 
the D2 line in Rb. However, the plasmon resonance is 
easily tuned by changing the geometry. In this work we 
consider the case of adding an inert core, such as SiC>2, 
to shift the plasmon resonance into the red [T3], as illus- 
trated in the inset to Fig. lc. Another approach would 
be to use nano-pyramids, where the dipole resonance is 
tunable by changing the aspect ratio of the pyramids [H] . 

Because this proposal requires the atom to be trapped 
within a few tens of nanometers from the sphere, there 
will be significant surface interactions. Van der Waals 
(vdw) forces can be overcome with modest laser power 
because of the strong plasmonic enhancement of the trap- 
ping laser near the sphere |15) . There are two dominant 
sources of heating and decoherence arising from incoher- 
ent transitions induced by the trapping laser or thermal 
magnetic field noise in the nanoparticle . The first effect 
scales as -7 Sl 2 /<5 2 , where 7 is the atomic linewidth, and 
is suppressed at large detuning. To estimate the effect of 
magnetic field noise we approximate the nanoshell as a 
current loop of radius and height a, thickness t, and resis- 
tivity p. Then the incoherent transition rate between hy- 
pcrfine states is ~ (gFPo^B) 2 kBT(aH/r 5 )/H 2 pr, where 
r is the distance of the atom to the sphere center, gp is 
the hyperfine g-factor, fig is the Bohr magneton, and T is 
the temperature |16| . Heating arises from both processes 
due to electromagnetic forces at the trap position which 



couple the motional states. The heating rate is reduced 
from the decoherence rate by a factor of 9h/2mr 2 ujT |17| . 

Figures led show the atomic trapping potential for a 
single sphere and an array, respectively. We numerically 
obtained the trapping potential in Fig. lc using Mie the- 
ory and the vdw potential was obtained using the meth- 
ods in Ref. jTS] . To solve for the trapping potential in the 
array in Fig. Id we approximated the scattered field from 
each nanoshell by a point dipole with a lorenztian polar- 
izability and solved for the total field self-consistently. 
Using the parameters in Fig. lc for trapping 87 Rb above 
a silver nanoshell at room temperature with £!o = 25 GHz 
(corresponding to ~ 10 8 J sat , where J sat ~ 1.7 mW/cm 2 ) 
and S = 25 THz, we estimate a trap depth of ~ 25 MHz 
and a trapping frequency of ~ 5 MHz. Both the mag- 
netic field noise and laser detuning limit the decoherence 
rate to ~ 10 Hz and the heating rate to ~ 1 Hz, meaning 
that the atom can be trapped for ~ 1 second. 

The controlled patterning of arrays of metallic 
nanoparticles can be done by lithographically defining 
such arrays in a top-down approach or through the 
controlled self-assembly of metallic nanoparticles, such 
as spheres, in a bottom- up approach 19-22]. In any 
nanofabricated system one must contend with disorder; 
the relevant disorder in this system occurs in the par- 
ticle positioning and particle formation. In lithographic 
approaches one can control the particle formation at the 
level of 1-2 nm [20] • In bottom- up, self-assembly ap- 
proaches it is possible to create large regions of well or- 
dered crystal with a finite density of point and line de- 
fects, much like a conventional solid [35] . Due to the local 
nature of the traps the disorder in the particle position- 
ing will not affect the trapping. Errors in the particle 
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formation can influence the trap by shifting the plasmon 
resonance and the field enhancement of each particle. To 
achieve consistent traps the fractional error in the plas- 
mon resonance should be smaller than its inverse qual- 
ity factor Q — lu sp /k, which for silver(gold) nanosphcres 
goes up to 80(20) depending the quality of the metal and 
finite size effects [12 US] . Currently, metallic nanoshclls 
can be made with a fractional error in the radius of less 
than 5%, which is comparable to the inverse of Q fM\ . 
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FIG. 2: Shows the scaling of the maximum tunneling in the 
lowest band, and the corresponding on-site interaction. Cal- 
culated using the Wannier functions for a sinusoidal poten- 
tial. (Inset) Energy dependent scattering length for two S7 Rb 
atoms on a single site as a function of the trap frequency. 



As a first example application of this system we con- 
sider a realization of the single-band Hubbard model in 
the novel regime of large atomic density pQ. Here the 
lattice spacing is limited only by fabrication technology 
and increasing surface interactions associated with scal- 
ing down the nanoparticle dimensions. As an example, 
Fig. Id shows that the vdw attraction can be overcome 
for a lattice spacing of 60 nm, which is within current 
fabrication limits. Figure 2a illustrates the scaling for 
the maximum tunneling in the lowest band and the cor- 
responding on-site interaction Uq, which scales as l/l 
in this limit due to the increased confinement [7J. In 
the supplementary material we show that the tunneling 
rate can also be tuned through appropriate polarization 
control [17] . 

In addition, these nanoscale traps reach a regime of 
atomic confinement where the ground state uncertainty 
becomes comparable to the free space scattering length 
(ao ~ 5 nm for 87 Rb). For two atoms in a 3D isotropic 
trap the two-body scattering problem can be solved ex- 
actly, leading to an effective scattering length a e g(ujT) 
which depends on the confinement energy [25, 26]. For 



example, a e g diverges when the scattering results in a 
it/2 phase shift on the wavefunction because the two- 
atom wavefunction will now coincide with the next ex- 
cited harmonic oscillator state. The inset of Figure 2 
shows the effective scattering length in the lowest vibra- 
tional state as a function of trap frequency [T7] . 

Disorder in the lattice has a strong effect on the Hub- 
bard model. The effect of disorder on the single-particle 
physics is well understood [57]; however, the interplay 
between interaction and disorder is an interesting new 
regime and will be considered elsewhere. 

As a second example, we now consider long range inter- 
actions within the plasmonic lattice, associated with the 
strong radiative coupling between the atoms and spheres 
[28] , In essence, this can be viewed as a strongly cou- 
pled cavity QED system. The dipolar coupling between 
the atoms and the near field of the sphere is given by 
g ~ nodo/eor 3 where do is the quantized dipole mo- 
ment of the sphere, which can be expressed in terms of 
the static polarizability a(0) as d = ,J huj sp a{0) / 2 [2"5] . 
Since the plasmons are overdamped the relevant coupling 
is given by the Purcell enhancement factor P — g 2 / 'ivy. 
The plasmon linewidth has contributions from radiative 
and ohmic losses. We approximate the ohmic damp- 
ing by its value in bulk, while the radiative damping 
rate is k^d^/iireoh ~ k 3 a 3 uj sp . Large spheres are ra- 
diatively broadened and, in this case, P ~ (fcr)~ 6 , while 
for small spheres P ~ Qa 3 /k 3 r 6 . In both limits, when 
r <C \/2ir ~ 100 nm the atoms enter the strong coupling 
regime P>1, see Figure [3^,. 

For a lattice of nanospheres, intersphere coupling is 
also present and leads to delocalized plasmon modes in 
the lattice [SHI US] • These modes can be used to medi- 
ate long-range atom-atom interactions. As an example, 
we calculate the interaction of two atoms through these 
modes in a ID chain of nanospheres. For each sphere 
in the chain we can write the self-consistent equation for 
their dipole moments as 

p n = a(u) (E n + N nm p m ) (3) 

where p n is the induced dipole moment of the nth 
nanopoarticle, E n is the incident field, and N nm is the 
3x3 matrix that gives the dipole field at site n due to the 
dipole at site m. In the ID system there are three orthog- 
onal bands corresponding to two sets of transverse modes 
where the dipoles are oriented perpendicular to the chain 
and one set of longitudinal modes where they are oriented 
parallel. Defining p q to be the qth eigenvector of N nm 
with eigenvalue D q , then the effective polarizability of 
the gth mode is a q l — a -1 — D q , i.e. p q = a q E q . For 
a Lorentzian polarizability the real part of D q gives the 
shift in the resonance frequency of the qth mode and the 
imaginary part gives the change in the linewidth. N nm 
is diagonalized by Fourier transform and if we neglect all 
but nearest neighbor terms D q = 27^ cosg — zfc 3 /67reo, 
where Nfa = Re(iYoi). 
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FIG. 3: a) Shows the cavity QED figure of merit g 2 /ivy 
with changing system size assuming the atom is trapped at 
a distance of twice the sphere radius. (Inset) Single atom 
trapped above a nanosphere acts as cavity QED system with 
atomic and cavity losses 7 and k, respectively, and a coherent 
coupling g. An optical field can be used to drive the cav- 
ity transition, b) Shows the fidelity for generating a ground 
state singlet state between two atoms on the lattice with their 
separation after optimizing the optical and microwave fields. 
The entanglement is generated through interaction with the 
collective plasmon modes, where we took a Q of 100 for the 
plasmon resonance. (Inset) The atom-plasmon array acts as 
a scalable cavity QED system. 



Let us consider atoms trapped above the ID array of 
spheres. The plasmonic modes can be adiabatically elim- 
inated using standard methods in quantum optics to give 
an effective interaction between the atoms with both co- 
herent and dissipative contributions .32] . For two-level 
atoms polarized parallel to the ID chain the atomic den- 
sity matrix evolution is 



a~}{ui a ) — 0. The first line in Eq. i] describes the coher- 
ent evolution and the second line describes the collective 
dissipation. Here we have neglected the contribution to 
the interaction from free-space radiative modes. 

Eq. 5 and 6 indicate that the coherent and dissipative 
contributions are equally strong when the atom and plas- 
mon are near resonant. Working far off resonance, how- 
ever, results in purely coherent dynamics, which can be 
used to implement long-range interacting spin models in- 
cluding frustration [33J GH] ■ Alternatively, the collective 
dissipative dynamics can be used to prepare correlated 
atomic states [35j . As an example, we now show how 
to directly prepare a ground state singlet between two 
atoms separated by large distances on the lattice. We 
take two ground states \g) and \s) and an excited state 
|e) which is coupled to \g) via an external field and only 
decays via the plasmons back to \g) (see inset to Fig. 3a). 
An external microwave field mixes the two ground states. 
To prepare the singlet state IS") = \gs) — \sg) we use a 
similar approach to Ref. [35] whereby the singlet state is 
engineered to be the steady state of a driven, dissipative 
evolution. We take a separation n such that cos q*n = 1, 
then in a rotating frame we can rewrite Eq. 4 

p = - l0n v[< + <4 e ] P - HnW] + v[4 e ]) P (7) 

where V[c]p = 1/2-jVc, p} - cpS and = 700 - 7n0 ~ 
700 (£ 3 /a 3 )n/Q for n < Q. The first term in Eq. 7 de- 
scribes the collective decay of the atoms into the plas- 
mons, while the second term is the residual linewidth due 
to plasmon propagation losses. In this case the dynamics 
can be mapped to a cavity QED system by identifying 
7on with the collective decay g 2 /k and S^ n with the free 
space decay 7. The two excited states \eg) and \ge) split 
into a superradiant state \eg) + \ge) and a subradiant 
state \eg) — \ge) with decay rates 27o„ + 8^ n and 5j n , 
respectively. 

The singlet preparation proceeds as follows. First, 
we selectively excite the subradiant transition \gg) to 
\ge) — \eg) by driving with a weak external laser field 
^ ~ Sjn <C 700, which we take to have a tt phase dif- 
ference on the two atoms. Second, in order to make the 
singlet state a unique steady state, we apply a global 
microwave field to mix the triplet ground states without 
affecting the singlet state. With this combination of ex- 
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where z is the position of the atoms above the sphere 
and q* = q* + i q* is the resonant wavevector such that 



while the triplet states are all excited at a much faster 
rate f2 2 /5j n . Additionally we assume that the state |e) 
can decay to |s) at a rate 6~/ n (which could be achieved 
e.g. by weakly mixing \e) with another excited state) [17] , 
In the resulting dynamics, the pumping rate into the sin- 
glet state is il 2 /5^ n , while the pumping rate back into the 
triplets is fL 2 /joo. The steady state of this process gives 
the singlet state with fidelity F = (S\ p \S) ~ 1 - 1/P' 
where P' = 7oo/^7n is the effective Purcell factor. Fig.j3|) 
shows the fidelity for two atoms with variable separa- 
tion as obtained from numerical simulation of the master 
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equation. Significant entanglement can be generated for 
atoms separated by many lattice sites. 

To measure the correlations in this system, an all op- 
tical approach could be realized by making the nanopar- 
ticle array in the near field of a solid immersion lens 
(SIL), which enhances the resolution beyond the diffrac- 
tion limit by a factor of n, the index of refraction of the 
SIL [37]. Combining a SIL with e.g. super resolution 
microscopy techniques would allow one to reach the req- 
uisite resolution of ~50 nm at optical wavelengths |38j . 

Our analysis shows that by combining cold atom tech- 
niques with nanoscale plasmonics, new regimes in con- 
trolling both the collective motion of atoms and atom- 
photon interactions can be reached. It may allow one to 
combine excellent quantum control of isolated atoms with 
nanoscale localization and strong interactions, thereby 
opening up exciting new possibilities for quantum con- 
trol of ultracold atoms. 
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